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We construct the most general form of axially symmetric SU (2)-Yang-Mills fields 
in Bianchi cosmologies. The dynamical evolution of axially symmetric YM fields in 
Bianchi I model is compared with the dynamical evolution of the electromagnetic field 
in Bianchi I and the fully isotropic YM field in Friedmann-Robertson- Walker cosmolo- 
gies. The stochastic properties of axially symmetric Bianchi I-Einstein- Yang-Mills 
systems are compared with those of axially symmetric YM fields in flat space. After 
numerical computation of Liapunov exponents in synchronous (cosmological) time, it 
is shown that the Bianchi I-EYM system has milder stochastic properties than the 
corresponding flat YM system. The Liapunov exponent is non-vanishing in conformal 
time. 
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1. Introduction 



The effects of anisotropy on the dynamics of the early universe have been a point of interest 
to cosmologists from time to time. This interest stems from the fact that by adding more 
degrees of freedom to any isotropic minisuperspace model one might hope to gain a better 
understanding of the behavior of the model generalized to the full superspace. Bianchi 
cosmologies with fluid sources are such models. The matter in these models is either a 
perfect fluid p4| or consists of massive or massless vector fields [25fl,||. 



There has also been interest in the study of homogenous source-free Yang-Mills fields as 
a dynamical system in the hope that a non-perturbative treatment might yield a better 
understanding of the vacuum state in YM theories, despite the fact that strong and weak 
interactions have no classical counterpart. The theory of these finite dimensional dynam- 
ical systems is dubbed Yang-Mills classical mechanics [ ^| . Similarly a non-perturbative 
mini-superspace Einstein- Yang-Mills (EYM) theory might eventually result in a better un- 
derstanding of the vacuum state of YM fields in the Planck regime. EYM cosmology is 
not new. There has been extensive work on various Friedmann-Robertson- Walker (FRW) 
cosmologies with a YM field source that has a stress-energy tensor of the form of a tracefree 
perfect fluid El , 1H , [0 , Q , [|0 . 



In this paper our aim is to relax the requirement of full isotropy. After adopting and refining 
a general scheme developed to construct YM fields on homogenous spaces, we examine, 
as a specific model, the dynamical properties of the EYM equations in axially symmetric 
Bianchi I cosmologies with an SU (2)-YM field. The organization is as follows: In section |2] 
after introducing the basic notation, we give a brief account of how invariant YM fields 
in Bianchi cosmologies with a given isometry group are constructed. This involves gauge 
fixing for both the space-time metric and the YM connection. The general field equations 
for invariant YM fields in Bianchi cosmologies are given in section ^. Then we use these 
equations to derive the evolution equations for axially symmetric YM fields in a Bianchi I 
cosmology followed by a brief review of how these equations are related to the known exact 
solution of axially symmetric electromagnetic fields in Bianchi I cosmologies and SU(2)-YM 
fields in FRW cosmologies. Section || contains a numerical analysis of the obtained EYM 
equations as a dynamical system, computation of Liapunov exponent and a comparison 
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with the flat space behavior. It is shown that surprisingly in synchronous time, the obtained 
EYM system has substantially milder stochastic properties than the corresponding flat space 
system. In conformal time, the Liapunov exponent is non-vanishing and the dynamical 
system is numerically less stable. 



2. Invariant YM fields in Bianchi cosmologies 

We consider Bianchi cosmologies where the space-time manifold is of the form R x E with 
a metric that admits an isometry group whose orbits are the space sections E t = {t} x E 
where E is a three-dimensional group manifold with a (t-dependent) invariant metric g. 
(This excludes the so-called Kantowski-Sachs solutions where E is not a group but only a 
homogeneous Riemannian manifold.) The space-time metric can then always be written in 



the form 22 



( g= -e ± (g ) e ± + g = -e 1 - ®e ± + g^t)? <g> e j , (i) 

where the 9 l (i = 1,2,3) are the components of the (left-invariant) Maurer-Cartan form on 
the group E and 9 ± = dt. If {e±,ei} is the (left-invariant) frame field dual to {^ _L ,6 |1 } then 

(4) 

the right-invariant vector fields are Killing vector fields of g (and g ) and they commute 
with the ej. 

The question of the most general form of the tensor fields on E invariant under certain group 
actions is extensively addressed in 0. Here we only briefly discuss a special case. We assume 
that the space-time admits a four- dimensional isometry group K whose orbits are the E t so 
that there is a one- dimensional isotropy subgroup K at each point. These so-called locally 
spatially isotropic cosmologies have been all been classified (see, for example, |TB|] ). The 
isotropy group K is then necessarily isomorphic to U(l) (as a one- dimensional subgroup of 
5*0(3)) and the metric g(t) can in all cases be chosen diagonal with two equal entries, 

{gab) = diag(6 x , 61,63), (2) 

say. In all cases but one (Bianchi III) it turns out that the action of U(l) on E is an 
automorphism of the group E that leaves the metric at the identity invariant from which 
it follows that K is a semidirect product of U(l) with E. The generator £^ of the isotropy 
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group at the identity then acts as an infinitesimal orthogonal transformation, and it follows 
that the commutation relations are 



= <?i£r (3) 

where explicitly 

= &,&]=£i, = o. (4) 

The question of invariance of the YM connection (or potential) in homogeneous spaces was 



addressed by Harnad et al. |12j . We give a short account of how such invariant connections 
are constructed in the case of Bianchi cosmologies. The complication is that an invariant 
YM connection is not necessarily constant in a left invariant frame (just as a Riemannian 
connection depending on the coordinate system does not necessarily vanish in Euclidean 
flat space), but any change in the field variables is merely due to a gauge transformation. 
However, it can be shown that for the YM potential 

( 4 ) 

A= A ± {x,t)8 x + A (5) 

in which A = Ai(x,t)8 l is the connection form on the homogeneous 3-space and A±(x,t) is 
a Lie algebra-valued scalar, there is always a gauge such that the Ai are only functions of 
time in a left invariant frame provided the 3-space is a group manifold (otherwise the 
9i are the pull back of Maurer-Cartan form components from the manifold of the isometry 
group to the homogenous space). This requires that all the Lie-algebra- valued fields that 
transform according to the adjoint representation (e.g. v4j_) have constant components in a 
left- invariant frame. Therefore the YM connection (|5|) reduces to 

{ A= A ± {t)6 L + AifyO 1 . (6) 

Several important facts should be mentioned regarding the YM connection constructed so 
far. 

(1) There is no local gauge freedom left in A+ and the remaining gauge freedom is global, 
i.e. only transformations of the form 

(Af)^^t)(Afh-\t), ( 7 6G) (7) 
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are allowed. (Here we have written Ai = Af Eg where {E^} is a basis of the Lie algebra q 
of the gauge group G. 

(2) With A±(t) 7^ 0, one can use the above global gauge freedom to make Af upper- 
triangular. The remaining six variables then represent the dynamical degrees of freedom of 
the SU(2)-YM field. 

(3) If an additional Killing vector field £^ generates an isotropy group U(l), it has a non- 
trivial action on the tangent space in view of the commutation relations The invari- 
ance of the YM connection requires the induced action of £^ on the cotangent space and on 
the Ai to be equivalent to a gauge transformation. 

To classify the possible .fT-invariant gauge fields systematically the following approach is 
needed |12|,[J15|. The equivalence classes of if-principal bundles P over E (where K is a 
Lie group that acts on P and acts via its projection by isometries on E) are in one-to-one 
correspondence with conjugacy classes of homomorphisms of the isotropy group Kq (= U(l) 
in our case) into the gauge group G (see These equivalence classes are well known from 

the investigations of spherically symmetric EYM-fields (p|,||17||,H) and are for K = U(l) K E 
and G = SU (2) classified by (nonnegative) integers n such that the n-th equivalence class is 
represented, for example, by 

A : 17(1) -> SU(2), e i(t> i-> e n ^ 3 (8) 

where tb = —icr B /2 G su(2) for B = 1,2,3 are used as a basis of su(2). On the other 
hand Wang's theorem (cf.|TI|) states that there is a one-to-one correspondence between the 
.fT-invariant G-connections on P and linear maps A : t — * g such that 

A o ad 2 = adx( z ) o A \fz e K . (9) 

and the connection components at the group identity can be chosen such that Ai = A(£j). 

Infinitesimally @ means in our case (K = U(l) tx E, K = U(l), G = SU(2)) that A*(^) = 
nr 3 and 

A([6,^])=n[A i ,r 3 ] (10) 
where we have put Aj = A(^) = A B i r B . Solving (]TD|), which becomes more explicitly 

A A r cl = ne A m A B l (11) 
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with cL as in (^) and (||), gives for n = 



for n = 1 



and for n > 1 





/o 





s\ 


A = 








e 




V° 





V 





/ 


a 


p 


o\ 


A = 




-a 


a 







V 








V 





(0 





o\ 


A = 













\o 





V 



(12) 



(13) 



(14) 



When these parameters a, (3, 7, 5, e are given as functions of t the YM connection is 
determined uniquely. 



3. Field equations for axially symmetric YM fields in Bianchi cos- 
mology 

We shall use units where 87r(Newton's constant) = (speed of light) = (YM coupling constant) 
1. The Yang-Mills field determined by A(t) in the gauge A± = is then 

f a = E^e 1 a e a + \B*e a a e\ (15) 

Et = At, (16) 
Bi = e A BC A^Af-A^. (17) 

where ' = C e± and the Lie algebra indices are raised and lowered with (5 a b) = diag(l, 1, 1) = 
{Sab) (the latter representing the invariant metric on su(2)). The YM equations are 

D^Ff^ = e A CD A Ci Ef - c r rs E As = 0, (YM constraint) (18) 
ZffjJ = Ef + K r r Ef - 2K{Ef + c\ s Bf r + \g ir J pq B A ™ + e A BC A Br B c ir = 0, (19) 



(3 



in which K ah = \ 'g a b and D„ is the 4-gauge-covariant derivative. Since 



T fw - F\ a F Av a - -g^F^F^ (20) 



the Einstein equations become 



(3) 

(K r r ) 2 — K rs K rs + R = E 2 + B 2 (scalar constraint), (21) 

K[c s rs + K r s c s ri = Bf r E r A , (momentum constraint) (22) 

Ky - 2K\K r3 + K r r K l3 + % = B A ir B A [ - EfE Aj + \g l3 {E 2 - B 2 ), (23) 
where E 2 = EfE r A and B 2 = \B*B r A s . 

For the metric (Q) the YM constraint is trivially satisfied if A has the form flT4[) . For the 
form (|12"1) it yields 

5e — 5e = £7 — £7 = ^7 — £7 = ^> (const. )e = (const. )5 = 7. (24) 
The YM constraint for (13) yields 



a/3 - a/3 = =^ a = (const.)/?. (25) 

The above equations show that after a time independent gauge transformation, ( |i~2D and 
(|I~3|) can be written in the following form, 

A = diag(0,0, 7 ), (26) 

A = diag(a,a,7). (27) 



Therefore modulo a gauge transformation, (^7|) is the most general form of an invariant 
SU (2)-YM connection in Bianchi cosmologies with a fourth Killing vector field obeying the 
commutation relations (|]). With this choice of the connection and the metric (fj) inserting 



c fc j = implies R\j = 0, i.e. the 3-space for Bianchi I models is flat. 



(3) 
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The evolution equations for axially symmetric YM fields in a Bianchi I cosmology (Bianchi I- 
EYM) are now 



a 2 , a 2 
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« + 7wT + «(f + -) = °> (28) 

(29) 

(30) 
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(32) 



bi 

in which (p8|) , (|29"D , (|3T]) and ( |32"D are the dynamical equations, ( pOf ) is the scalar constraint 
and ' = (where is t is the synchronous time). 

We consider first two special cases. 

Electro-magnetism: With a = the case (|i~3"D reduces to ([Mf) . With (d/dt) = (y/bs/bi)(d/dr), 
the general solution to the YM equations is 7 = c\T. Subtracting (|30|) from (|3T|) and adding 



(j30[) to (|32"D , respectively, with a time reparametrization (d/dt) = y/bs(d/dT') gives the 
solution 

b t = (cqt'+v 7 ^) 2 , (33) 

6 3 = ^_ (34) 

c r' + v^ (c r' + v ^) 2 

in which cq,Ci,B and A are the integration constants. This solution is equivalent to the 
known solution of the Einstein-Maxwell equations for an electromagnetic field in an axially 



symmetric Bianchi I universe . The energy- momentum tensor in an orthonormal frame is 
(Tjf) = diag(— p, p, p, —p) in which p is the matter energy density. Heuristically, the positive 
principal pressures in directions 1-2 and negative pressure in direction 3 explain why such 
a universe evolves as equations fl3"3"|)and fl3"lp indicate. During any expansion in direction 3 
energy is transferred from the gravitational field to the EM field whereas in any expansion 
in the 1-2 directions, energy is transferred from the EM field to the gravitational field. 
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However there is no potential energy associated with the gravitational field. Therefore there 
is an expansion in 1-2 directions and any expansion in direction 3 can not be sustained for 
a long time. In this model the Ricci tensor uniquely determines the EM field tensor up to a 
constant duality transformation. 

Isotropic case: Imposing spherical symmetry such that Kq = SU(2) requires a = 7 and 
b\ = 63 in which case the EYM equations reduce to those for a SU(2)-YM field in a FRW 
cosmology. In conformal time the EYM equations are given in ||10|. The solution for the 
YM field variables is given by elliptic integrals. The energy-momentum tensor is that of a 
radiation perfect fluid with energy-momentum tensor (T^) = diag(— p, p/3, p/3, p/3) and the 
geometry is that of a Tolman universe in which the space-like hypersurfaces of homogeneity 
are flat. In synchronous time b\ = b 3 = C\t + C2 where C\ and C2 are integration constants. 
In this particular example, one can easily show that any axially symmetric YM connection 
must necessarily be spherically symmetric. A comprehensive treatment of Emstem-S77(n)- 
YM system in FRW cosmologies is in preparation. 

To further facilitate the analysis of the dynamical system (p8|) -(p2|) one can use the Hamil- 
tonian 



«'i-(j + «V)]/A 



(35) 



in which Pi,P 3 ,P a and P 7 are the momenta conjugate to bi, b 3 , a and 7 respectively. Now 
the system (f2~§|)-(|32|) can be written in the equivalent form 
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(36) 



and the constraint H = 0. One can convert the above system of equations into polynomial 
form by the time reparametrization dr = dt/(bi\/bs) and a transformation s» = Pib^ i = 
1,3. It turns out that the only set of equilibrium points of this system in the physical 
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region b\ > 0, 63 > is the invariant submanifold si = S3 = a = P a = P 7 = which 
corresponds to flat space. Correspondingly the dynamical equations of motion in conformal 
time drj := {b\bz)~ 1 ^ % dt are derived from the Hamiltonian -^conformal = {hi/b^Y 1 /^ y/bsH. 
The above system is invariant under the group of scale transformations a — > col, 7 — ► 
cry, P a — > c 2 P a ,P 1 — > c 2 P 1 ,Pi — > cPi,P 3 — > cP%,bi — > c 2 b\,bz — > c 2 ^. One can use this 
symmetry to reduce the order of the above system from eight to six. However, due to the sin- 
gular nature of the transformation, the resulting system is not suitable for numerical analysis. 
The energy-momentum tensor in an orthonormal frame is (Tjf) = diag(— A — B, B, B, A — B) 
where 

A ~ hb, ' (3?) 



2b\ 



B = (38) 



Contrary to the electromagnetic and fully isotropic cases the principal pressure in direction 3 
does not have a definite sign. As will be seen from numerical investigations, the numerators 
of both expressions (^)and fl38|) have the same order of magnitude. However, any decrease 
in 63 will cause the positive term in T| to dominate and (note the discussion after (|34])) 
63 starts to increase. Hence, generally speaking, one would expect both b\ and 63 to be 
increasing functions of time. In fact, we used a fifth order Runge-Kutta integrator to 
integrate the system (|36D and computed the Hamiltonian constraint to check the accuracy 
of the numerical integration. 

As figure 2 and the YM equations indicate, the general behavior of the YM fields is that of 
two coupled anharmonic oscillators with time- dependent frequencies. The behavior of the 
YM field variables in the above system resembles the dynamical properties of homogeneous 
YM fields in flat space known as Yang-Mills Classical Mechanics (YMCM) @. 



4. Axially symmetric YM fields in flat space and regularizing ef- 
fects of gravitational self-interaction 

A full analysis of the dynamical system ( |36D is an insurmountable task. Therefore we de- 
cided to start our analysis from the simpler system of axially symmetric YM fields in flat 
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Figure 1: The behavior of metric variables for initial conditions b\ = 0.50,63 = 11.12, P3 = 
0.03, a = 7 = 0.1, P a = -0.67, P 7 = 0.03. The inserted figure covers < t < 100 and 
indicates an initial decrease in 63 which quickly reverses direction. 



space. Fortunately, the procedure described in section |2| encompasses the gauge fixing for 
homogeneous YM fields in flat space. It is well known that YMCM has stochastic properties 
P^|,[23|,|j6|,[14|. In these models the reduction from the full space of dynamical variables 



to lower dimensions to make the dynamical evolution tractable is via some ansatz. In our 
model, the reduction is an inevitable consequence of the space-time symmetry. 

The two dimensional flat system A\ 7^ 0, A2 7^ and the three dimensional flat system A\ 7^ 
0, A2 7^ 0, A3 7^ 0, all other components vanishing, have been extensively covered in |23| and 
|J. The stochastic character of these systems is demonstrated by numerically computing the 
Liapunov index. Such numerical computations are achieved by the simultaneous integration 
of the first order system x = V(x) and the linearized first order system w = M(x).w in 
which w is the perturbation vector connecting two nearby trajectories and M(x) is the 



Jacobian matrix of V(x) [fL8||. The Liapunov exponent is defined as 



1, w(t) , s 

a = lim - In ^Hr- 39 

t^oct |w(0)| V ^ 

We follow the procedure explained in to compute the Liapunov index for axially symmetric 
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Figure 2: The behavior of YM field variables in axially symmetric Bianchi I-EYM cosmology 
for the initial conditions b\ — 2, 63 = 1, P3 = a — 7 = P a — 0.1, P 7 = 0.2. The oscillations of 
7 are characterized by higher amplitudes and lower frequencies. 

YM fields first in flat space and later in a Bianchi I cosmology. 

The dynamics of axially symmetric YM fields in flat space is governed by the system 



which correspond to a set of two strongly coupled oscillators with varying frequencies and 
amplitudes. These equations are derived from the Hamiltonian 



The system describes the motion of a point particle moving in a potential well U = a 2r y 2 + 
a 4 /2 with two open channels in the directions of positive and negative 7 (figure §). In these 
channels the term quartic in a in U is much smaller than the term quadratic in a. Therefore 
the behavior of the point particle in each channel is basically the same as that of the point 
particle in the two dimensional Hamiltonian system 



a + a(7 2 + a 2 ) 
7 + 2a 2 7 



0. 



(40) 
(41) 




(42) 



tt 1 / • 2 • 2 1 2 2\ 

H= -{a 7 +a 7 ) 



(43) 
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treated in JHJ. The potential barrier in this system has open channels both in a and 7 
directions and such systems have been extensively studied because of their relation to the 
plasma confinement problems. Unless a = or 7 = 0, as the particle moves deeper and 
deeper into the 7 channel, say, the frequency of oscillations in a increases while the amplitude 
decreases. However, at a finite value of 7, 7 = at which point the particle returns to the 
a ~ 7 region. In this system, stochastic regions occupy a significant portion of the phase 
space and the regular region is limited to a very small (~ 0.005%) region of the phase space 

H- 

7 



a 



Figure 3: The behavior of a particle in potential well U = a 2r y 2 + a 4 /2. 

The behavior of a particle in a system with potential barrier U is basically the same. However, 
because of the lack of the existing channels in directions a, oscillations of 7 are characterized 
by larger amplitudes and smaller frequencies. 

Following we numerically computed the Liapunov index for the system (pE0"|) and (pUP 
for randomly selected initial conditions satisfying (f42|) with H — 1. It turns out that the 
Liapunov index for this system is positive and is of the same order of magnitude as the 
one for the system (pE3"D (see figure |3j). Numerical investigations for randomly selected initial 
conditions indicate that in this system, stochastic regions occupy a large portion of the phase 
space also. 
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Figure 4: In of the Liapunov exponent vs. In of the evolution parameter for axially symmetric 
YM fields in flat space and in Bianchi I-EYM cosmology for initial conditions as in figure 
2. This was a typical behavior for randomly selected initial conditions, (s) and (c) refer to 
synchronous and conformal time respectively. 



There are at least two problems associated with generalizing the study of the stochastic 
properties of the flat space model to axially symmetric Bianchi I-EYM model represented by 
the system (|36l). One is related to the strongly coupled nature of the ODE system and the 
higher number of degrees of freedom which are known to cause sophisticated stochastic phase 
space properties like Arnold diffusion fL8| . Numerical investigations (see figure [D point to a 
non-compact phase space. Thus we can say that axially symmetric Bianchi I-EYM systems 
are not globally ergodic. However, we do not rule out the existence of ergodic components. 

The other problem is related to the inherent gauge dependence in the definition of Liapunov 
exponent and the non-existence of a satisfactory gauge covariant definition of chaos in general 
relativity. It is known that in Mixmaster models the positivity of Liapunov exponents 
depends on the choice of time reparametrization However, in Mixmaster cosmology, the 
stochasticity is associated with the behavior of the metric variables in the vicinity of the 
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cosmological singularity where cosmological time is not well defined. In Bianchi I-EYM any 
ergodicity, if there is any, is mainly in YM field variables, far away from the cosmological 
singularity 

Similar to the flat space scenario, we calculated the Liapunov exponents in both synchronous 
and conformal time of Bianchi I-EYM system for randomly selected initial values (see figure 
H). The vanishing of Liapunov exponents in synchronous time point to a dynamical system 
in which the stochastic regions, if there are any, occupy a much smaller portion of the phase 
space. However , it also underlines the known fact that the Liapunov exponent is sensitive 
to time reparametrization. The Liapunov exponent is non-vanishing in conformal time with 
a value smaller than the corresponding flat space model. At this point we would like to add 
that it is more difficult to preserve the constraint in the conformal time and the numerical 
stability of the dynamical system is substantially enhanced in the synchronous time. 

Following H we computed the correlation between the initial and final values of YM field 
variables as an indication of a particular statistical independence in the dynamical evolution 
of the field variables. As figure |5] demonstrates, after a large time evolution, there is a loss 
of correlation between ^ and 7/ (respectively and 07). Therefore 7* (a*) and 7/(07) can 
be regarded as two stochastically independent random variables. 
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Figure 5: The correlation between 7^ and ji in terms of time for the initial conditions 
b x = b 3 = 1, P x = 2, P 3 = 1, P a = -.28. In the left figure a = 0.1, 0.15 < 7 < 17.25 and in 
the right figure 7 = 0.1, 0.15 < a < 1.56. The large time behavior indicates that there is no 
clear correlation between 7/(07) and 7i(oj). 




8 -0.5 - 



15 



5. Conclusion 

We systematically derived the most general form of the YM connection and the EYM field 
equations in Bianchi cosmologies with a four-dimensional isometry group in which the Killing 
vector fields obey the commutation relations (f|). For the simplest of Bianchi cosmologies, 
namely Bianchi I, we investigated the resulting dynamical system. In doing so, one realizes 
that there is little hope to find an exact solution. Numerical integration suggests a non- 
compact phase space and oscillatory behavior for the YM field variables. However, one can 
easily use the above mentioned scheme to construct invariant YM connections in flat space. 
There has been extensive work on the dynamical properties of homogeneous YM fields in flat 
space (YMCM) which are known to have stochastic properties. We used some methods to 
investigate chaos in YMCM (i.e numerical computation of Liapunov exponent) to see how 
gravitational self-interaction can affect the stochastic behavior. It turned out that the system 
with gravitational self-interaction has milder stochastic properties. We hope to extend this 
work to other Bianchi cosmologies. 
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